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A  STUDY  OF  THE  FEASIBILITY  OF  INCREASING  CIRCULATION 
IN  THE  WATERWAYS  OF  THE  COLLIER  BAY  AREA 
MARCO  ISLAND,  FLORIDA 


.Purpose  of  the  Study 

This  study  was  designed  to  determine  the  feasibility  of  improving 
the  circulation  of  the  Landmark  and  Castaways  Waterway  systems  at  Marco 
Island.  In  view  of  the  fact  that  many  of  the  waterways  had  already  been 
built,  only  minor  changes  in  configuration  could  be  considered. 

Description  of  the  Area 

Marco  Island  is  located  on  the  Gulf  Coast  of  Florida  at  Latitude 
25°56'N.  The  island  itself  is  approximately  four  miles  long  and  is 
isolated  from  the  mainland  by  the  Marco  River.  Major  inlets  consisting 
of  Big  Marco  Pass  to  the  north  and  Caxambas  Pass  to  the  south  provide 
good  exchange  between  the  River  and  the  Gulf.  The  tide  is  well  mixed 
and  has  a  mean  range  of  1.7  feet  and  a  spring  range  of  2.6  feet. 
Description  of  the  Systems 

The  waterways  systems  are  shown  on  Figure  1.  The  Landmark  system  has 
a  planned  surface  area  of  approximately  2200  acres  and  is  tributary  to 
Caxambas  Pass  on  the  south.  The  longest  distance  to  Caxambas  Pass  is 
11,000  feet. 

The  Castaways  system  has  an  area  of  2400  acres  and  drains  into  the 
Marco  River  near  Big  Marco  Pass  via  Collier  Bay.  Its  maximum  distance  to 
the  Marco  River  is  9,500  feet. 

The  landward  ends  of  the  two  waterways  are  separated  only  by  State 
Road  92  and  a  tier  of  lots  on  both  sides,  a  total  of  350  feet. 


Approach  to  the  Problem 


In  attacking  this  problem,  it  was  immediately  evident  that  the 
proximity  of  the  landward  ends  of  the  waterways  might  be  used  to  ad¬ 
vantage  by  connecting  them  at  one  or  more  points.  A  quick  study  of 
the  records  of  tide  gauges  located  near  the  mouths  of  the  system 
showed  a  phase  difference.  This  indicated  that  a  differential  head 
would  exist  across  State  Road  92.  This  head  would  be  available  to 
promote  interchange  between  the  two  systems.  The  actual  extent  of 
the  interchange  could  be  predicted  by  means  of  a  numerical  model. 

The  Numerical  Model 

This  model  was  based  on  the  application  of  Manning's  formula  to 
a  series  of  separate  reaches  in  the  system.  Reaches  are  selected  so 
that  their  geometry  can  be  considered  to  be  uniform  throughout.  For 
each  reach,  a  storage  area  and  a  resistance  coefficient  is  determined. 

The  following  relationships  will  then  hold  for  the  nth  reach  for 
continuity  and  energy  loss : 


where : 

3 

Qn  ■  flow  into  the  reach  (ft  /s) 

3 

Qn+1  ■  flow  out  of  the  reach  (ft  /s) 

2 

Ar  "  storage  area  (ft  ) 

All  ■  increase  in  surface  elevation  (ft) 
At  -  interval  of  time  (s) 


The  heed  loss  in  the  reach  is  related  to  the  flow  by 


T1  "  -  k  /  Q  +  Q 

'n  'irt-l  n  /  _n _ 2 


where: 


T|  -  water  surface  elevation  at  entrance  (ft) 
n 

-  water  surface  elevation  at  exit  (ft) 

2 

k  *»  resistance  coefficient  for  the  reach  (s/ft  ) 


where:  k  is  derived  from  Manning  s  formula 

u  -  1,486  R2/3  Sl/2 

n  — 2 -  n  n 

n 


Vs  ^n  +  Vl)  -  »>n  an 


where :  _  3 

Q  *  the  root  mean  square  discharge  in  the  reach(ft  /s) 
n 

u  »  velocity  averaged  vertically  and  horizontally  (ft/s) 
n 

1.486  -  a  conversion  factor  from  metric  to  English  units. 

(3.281  ft/m)l/3 

n  -  Manning's  roughness  coef f icient(ft) 

R  -  the  hydraulic  radius  (ft)  which  is  the  cross  sectional 
n  area  of  the  channel,  a  (ft2)divided  by  the  wetted  perimeter, 
Pn(ft) 

S  »  the  hydraulic  gradient(f t/f t) 
which  for  uniform  flow 

—  ATJ  4-  Ln  where  Lr  -  the  length  of  the  reach(ft) 


s 


Substituting  in  Equation  (4)  and  squaring  we  have 

10/3 
”4/3 


Q2  =  (  1.486  )2  an 
n  V  ~TT"  '  “4 


£1 

L_ 


which  gives: 


2  P*'3  t, 

n _ n _  n 

2,208  a  10/3 


(5) 


(6) 


Obviously  in  an  area  affected  by  tide  kn  will  vary  with  time  as  pn 
and  a^  vary  with  the  water  surface  elevation.  In  the  present  model,  this 
variation  was  found  to  be  insignificant  due  to  the  large  depth  (20  ft)  of 
the  canals;  so  it  was  ignored. 

Figure  2  shows  a  schematic  diagram  of  the  waterways  systems  as 
modelled. 

The  boundary  conditions  are  prescribed  by  the  water  surface  elevations 
at  the  mouth  and  the  geometry. 

In  analyzing  these  elevations,  the  raw  data  from  tide  gauges  is  smoothed 
to  eliminate  transient  effects  such  as  wind  set  up  by  utilizing  a  tidal 
component  analysis  involving  a  least  squares  fit  to  isolate  the  mean  tide 
level  and  sine  and  cosine  coefficients  for  the  four  principal  frequencies 
of  tide  producing  forces  according  to  the  formula: 

4 

Tl  ■  Tl  +  £,  (S  sin  o  t  +  C  cos  a  t) 

'  o  n-1  n  n  n  n 


m  cm  oo 


O  oo  o\ 


SCHEMATIC  DIAGRAM  OF  NUMERICAL  MODEL 


where: 


T]  =  tide  level  at  time  t 

11  =  mean  tide  level 

‘o 

■  coefficient  for  the  sine  term  for  the  nth  frequency 

Cn  =  coefficient  for  the  cosine  term  for  the  nth  frequency 

0  =  frequency  or  2tt/T  where  x  is  the  period  of  the  nth 

n  ,  n  n 

frequency 

For  this  study  the  analysis  included  components  having  periods  as 
shown:  M2-12.42  hours,  S2-12.00  hours,  Kl-23.93  hours  and  01-25.82  hours. 

The  components  are  used  to  compute  typical  spring  and  neap  tides  to 
be  used  in  the  model.  In  the  case  of  this  study  hourly  values  were  com¬ 
puted  for  26  hour  periods  at  spring  and  neap  tide. 

In  operating  the  model,  the  first  approximation  for  discharges  is 
based  on  a  frictionless  condition.  The  head  loss  and  storage  volume 
is  then  recomputed  by  Equations  (1)  and  (2)  on  the  basis  of  these  flows 
to  determine  a  second  approximation.  This  is  averaged  with  the  preceeding 
value  and  the  process  is  repeated  until  a  reasonable  tolerance  is  obtained. 
Application  of  the  Model 

On  the  basis  of  using  the  tidal  differences  in  the  two  waterway  systems 
to  increase  flushing,  the  model  as  shown  on  Figure  2  was  revised  to  pro¬ 
vide  a  channel  20  feet  wide  and  10  feet  deep  connecting  Reach  No.  105 
(Castaways  Waterway)  of  the  Castaways  system  with  Reach  No.  510  (Shenandoah 
Waterway)  of  the  Landmark  system.  A  similar  channel  was  provided  to 
connect  Reach  No.  107  (Copperfield  Waterway)  with  Reach  No.  508  (Landmark 
Waterway).  These  connections  are  shown  by  dotted  lines  on  Figure  2. 

These  channels  were  found  to  have  a  resistance  coefficient,  k  -  100,000x10”* 


j 
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considerably  more  than  the  total  resistance  of  all  channels  in  the 
network. 

Operation  of  the  revised  model  utilizing  existing  data  from  the 
David  Key  and  Yacht  Club  tide  gauges  and  allowing  25%  of  the  velocity 
head  for  entrance  and  exit  losses  showed  that  the  discharge  through 
the  canals  would  be  more  than  trebled  on  both  spring  and  neap  tides. 

The  situation  was  further  improved  by  simulating  gates  at  the 
channels  to  allow  flow  only  to  the  north.  Under  these  conditions,  it 
was  shown  that  the  net  rate  of  exchange  per  24  hour  period  would  exceed 

130  x  10^  cubic  feet  on  a  spring  tide  and  90  x  10^  cubic  feet  on  a  neap 
tide. 

Since  the  tidal  data  used  could  not  be  assumed  to  reflect  the 

exact  conditions  in  the  waterways,  two  recording  gauges  were  installed 

in  the  waterways  themselves  and  operated  during  the  period  1-26  April 

1973.  These  gauges  were  accurately  referenced  to  each  other  and  to 

Sea  Level  Datum,  1929,  by  differential  leveling. 

These  data  were  then  reduced  and  tides  were  reconstructed  for  10 

April  1973  to  produce  a  neap  tide  and  for  17  April  to  produce  a  spring 

tide.  Plots  of  these  are  shown  on  Figures  3  and  4. 

The  model  was  then  operated  utilizing  the  reconstructed  tides. 

Conditions  with  and  without  gates  were  assumed  and  volumes  of  net  flow 

toward  the  north  were  computed.  From  these  volumes  and  the  total 

6  3 

volume  of  the  waterway  systems:  Castaways  -  145  x  10  ft  ,  Landmark  - 
6  3 

131  x  10  ft  ,  the  exchange  time  for  the  total  volumes  was  computed. 
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FIGURE  3  -  TIDAL  RELATIONS  FOR  LANDMARK  AND  CASTAWAYS  WATERWAYS 


FOR  NEAP  TIDE 
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FIGURE  4  -  TIDAL  RELATION  FOR  LANDMARK  AND  CASTAWAYS  WATERWAYS 

FOR  SPRING  TIDE 


The  results  of  operating  the  model  are  shown  in  detail  in  the  Appendix. 
A  summary  is  shown  on  Table  1. 


1.  Without  Gates 


TABLE  1 


Net  North  Flow 

3  -6 

(ft  /day  x  10  ) 


Time  to  Exchange 
Total  Waterway 
Volume 

_ (days) _ 


a.  Castaways  System 

(1)  Neap  Tide 

81 

1.80 

(2)  Spring  Tide 

33 

4.39 

b.  Landmark  System 

(1)  Neap  Tide 

87 

1.51 

(2)  Spring  Tide 

21 

4.49 

2.  With  Gates 

a .  Castaways  System 


(1)  Neap  Tide 

112 

1.31 

(2)  Spring  Tide 

96 

1.53 

Landmark  System  * 

(1)  Neap  Tide 

117 

1.12 

(2)  Spring  Tide 

92 

1.43 

Conclusions 

As  a  result  of  this  study,  we  conclude  that  the  circulation  of  the 
Castaways  and  Landmark  Waterway  systems  can  be  greatly  improved  by  the 
construction  of  two  concrete  channels  20  ft.  wide  by  10  ft.  deep  by  350 
feet  long  under  State  Road  92.  These  channels  will  provide  total  ex¬ 
change  in  one  to  three  days  which  will  provide  exceptionally  good 
water  quality  as  compared  to  the  configuration  as  currently  planned. 

The  installation  of  gates  that  will  alow  flow  only  to  the  north 


will  further  improve  the  circulation,  but  the  cost  of  constructing 


maintaining  and  operating  these  gates  will  outweigh  the  benefits 
gained . 

Recommendations 

We  recommend  as  follows : 

1.  That  the  westernmost  channel  connecting  the  existing  Landmark 
and  Copperfield  Waterways  be  constructed  as  soon  as  practicable. 

2.  That  after  this  is  constructed  the  flows  in  the  existing  system 
be  analyzed  in  detail  to  check  the  authenticity  of  the  model. 

3.  That  a  channel  connecting  Shenandoah  Waterway  with  Castaways 
Waterway  be  included  in  future  plans  for  developing  the  Collier  Bay 
Tract. 

4.  That  the  design  of  both  channels  be  modified  on  the  basis  of 


results  of  the  analysis  recommended  in  Paragraph  2  above.  At  this  stage 
it  may  be  considered  desirable  to  add  tide  gates. 
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APPENDIX 


RESULTS  OF  MODEL.  RUNS 
CONTENTS 


Explanation  of  Format 

I 

TABLE  A  -  Castaways  Waterway 

Without  Tide  Gate 

II 

TABLE  B  -  Landmark  Waterway 

Without  Tide  Gate 

III 

TABLE  C  -  Castaways  Waterway 

With  Tide  Gate 

IV 

TABLE  D  -  Landmark  Waterway 

With  Tide  Gate 

V 

14 


Explanation  of  Format 

These  tables  show  flow  conditions  for  neap  and  spring  tides, 
columns  are  as  follows: 

Column 

1  &  7  Time  in  hours  after  0000  hrs .  on  1  April  1973 

2  &  8  Tide  height  above  1EL  1929  datum  at  mouth  of 

sys  tem 

3 

3  &  9  Flow  in  ft  /s  through  the  channel.  Flow  to  the 

north  is  positive. 

4  &  10  Flow  past  Station  1  marked  C  1  and  L  1  on  Figure 

2  for  Castaways  and  Landmark  respectively 

5  &  11  Flow  past  Station  2 

6  &  12  Flow  past  .Station  3 
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TABLE  D  -  LANDMARK  WATERWAY  WITH  TIDE  GATE 


